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Abstract 

This is a review of selected work on the one- and multi-dimensional 
random Burgers equation (burgulence) with emphasis on questions gen- 
erally asked for incompressible Navier-Stokes turbulence, such as the law 
of decay of the energy and the pdf of velocity gradients. Most of the ma- 
terial is devoted to decaying (unforced) burgulence. For more details see 
the Table of Contents. 
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1 Introduction 



These lectures are about the d-dimensional Burgers equation 

d t v + (v • V) v = iA7 2 v, v = -V-0. (1) 

Note that the constraint that v be derived from a (velocity) potential if> is 
trivially satisfied if d = 1. The word "burgulence" , as we use it here, is a 
contraction of "Burgers" and "turbulence" . It means "the study of random 
solutions to the Burgers equation" . The randomness may arise because random 
initial conditions vo = — V^o are given or because a random driving force 
f = -VF is added to the r.h.s. of (0), or both. When f = one speaks about 
"decaying burgulence" . 

In the thirties when the Dutch scientist J.M. Burgers introduced the equation 
in the one-dimensional case, he hoped to contribute to the study of turbulence 
with a simple model which, obviously, has a lot in common with the Navier- 
Stokes equation: 

• same type of advective nonlinearity 

• presence of a diffusion term from which a Reynolds number may be defined 

• many invariance and conservation laws in common: invariance under 
translations in space and time, parity invariance, conservation of momen- 
tum and energy (only for v = and d= 1). 

Such hopes appeared to be shattered when, in the fifties, Hopf Q and Cole || 
discovered - some say rediscovered - that the Burgers equation can actually be 
integrated explicitly (we shall return to this matter later). Indeed, an important 
property of the Navier-Stokes equation, not shared by the Burgers equation, 
is the sensitivity to small changes in the initial conditions in the presence of 
boundaries or driving forces and at sufficiently high Reynolds numbers. Hence, 
the Burgers equation is not a good model for one of the most important aspects 
of turbulence: the spontaneous arise of randomness by chaotic dynamics. 

In spite of this there has been a strong renewal of interest in the Burgers 
equation, starting in the eighties, for a variety of reasons which we shall now 
explain. As a quantitative measure of the current interest, Table [j] shows some 
web-based statistical figures on the number of hits as of August 2000 (Google is 
an all-purpose search engine and "Los Alamos" stands for the nlin (ex-chao-dyn) 
preprint archive) : The Burgers equation, which obviously describes a compress- 





Navier-Stokes equation 


Burgers equation 


Google 


15000 


4000 


Los Alamos 


100 


75 



Table 1: Web-based statistical data 



ible flow (in one dimension there exist only trivial incompressible flows), has 
found many applications in nonlinear acoustics and other nonlinear wave prob- 
lems. A review may be found in Ref. H. 
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1.1 The Burgers equation in cosmology 

The Burgers equation has found interesting applications in cosmology, where it 
is known, in one instance, as the "Zel'dovich approximation" |jj and, in another 
instance, as the "adhesion model" Q. Here, we shall give a brief introduction 
to how the Burgers equation arises in cosmology. More details may be found in 
Refs. H 0, ||, Just after the baryon-photon decoupling in the early Universe, 
there may have been a rarehed medium formed by collisionless dustlike particles 
without pressure, interacting only via Newtonian gravity ||. The gravitational 
potential is then determined from the fluctuations in mass density by a Poisson 
equation. Limiting ourselves to the case of a single type of matter, we can 
schematically write the acceleration of a fluid particle as follows: 

acceleration = pressure + viscous + expansion + gravit. 

term term term term 

ocv in 

d t v + v • Vv negligible ? comov. coord. 

On the left hand side (l.h.s.) we recognize the familiar terms of the Burgers 
equation. The pressure is usually neglected because the matter is very cold. We 
shall come back to the viscous term later. The expansion term, proportional 
to the velocity, arises because the equation is written in a frame comoving with 
the expansion of the Universe. 

It turns out that when the problem of self-gravitating gas in an expanding 
universe is examined in the linear approximation (small density fluctuations) an 
instability is obtained in which the dominant mode has the following properties 
§|: 

• it is potential (v = — Vi/'), 

• the expansion and gravitational terms cancel. 

In 1970, Zel'dovich Q proposed to extend these properties into the nonlinear 
regime where density fluctuations become strong and mass condensation forming 
large-scale structures appear. Furthermore, this "Zel'dovich approximation" is 
exact in one dimension, irrespective of the strength of fluctuations. Clearly, 
in the Zel'dovich approximation each fluid particle is just moving in a straight 
line with constant velocity (after a suitable nonlinear change of variable of the 
time). Just like a family of straight light rays forms generally caustics along 
which the intensity is infinite, the material particle lines form singular objects 
along which the mass density is infinite. Arnold, Shandarin and Zel'dovich 
[ fl0| studied the various kind of singularities which can form in this way, to 
see if they could account for known large-scale structures such as galaxies and 
clusters. Observations and numerical simulations have now revealed that large- 
scale structures are much simpler than the mathematical objects generated in 
a caustic-type theory. 

Consider for example Fig. |l| which shows a thin slice of a simulated uni- 
verse using the so-called ACDM model (cold dark matter with cosmological 
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Figure 1: iV-body simulation by the Virgo Consortium (see 
http://star- www. dur.ac.uk/^frazerp/virgo/virgo. html). The simulation 

has 256 3 particles and was done on two large Cray T3D parallel supercom- 
puters at the computing centers Garching (D) and Edinburgh (GB). The 
brightness is proportional to the log of the density of the particles. 
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Figure 2: Mass density field from a two-dimensional 512 2 simulation of the 
Burgers equation with random scale- free initial data (from Ref. 
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constant). The line- like and node-like features on this slice are actually sheets 
and filaments embedded in the three-dimensional space. Sheets (walls), fila- 
ments and nodes (clusters) are the most common structures observed in such 
simulations. As was shown by Gurbatov and Saichev || these are precisely 
the structures obtained if one modifies the Zel'dovich approximation by requir- 
ing that particles should not cross but rather adhere. This adhesion model is 
just the three-dimensional Burgers equation (|l]), taken in the limit of vanishing 
viscosity. Numerical experiments indicate that the adhesion model reproduces 
quite well the early skeleton of large-scale structures in N-body numerical sim- 
ulations (see, for example, Figs. 6a and 6b of Ref. (nj]). Since dark matter is 
essentially collisionless, it is not clear at the moment what is the physics behind 
this agreement which seems to require some viscosity-generating mechanism to 
prevent - or dramatically slow down - particle crossing. Furthermore, the adhe- 
sion model cannot cope with many important aspects of gravitational dynamics. 
For example, in .ZV-body simulations, one frequently observes the collapse of a 
filament into an isolated node (cluster). As we shall see, there is nothing of this 
sort in Burgers dynamics. 

1.2 The Burgers equation in condensed matter and sta- 
tistical physics 

The Burgers equation arises in a number of condensed matter and statistical 
physics problems and even in non-physics problems such as vehicular traffic (for 
review see Ref. 0[) A frequently studied problem is the Kardar-Parisi-Zhang 
or KPZ equation [jj (see also Ref. @) 

dttp = \ |W>| 2 + i/V 2 V + F, (2) 

which appears in studying the motion of an interface under deposition. Here, 
ip is the vertical displacement of the interface as a function of d — 1 horizontal 
coordinates and of the time. It is immediately checked, by taking the horizontal 
gradient of (||), that one obtains the Burgers equation ([!]) with an additional 
forcing term f = — VF. Burgers equation also arises in studying directed poly- 
mers (see, e.g., Refs. jl5[ |l6|), but with the time variable now interpreted as a 
space variable in the direction of main extension of the polymers. On all these 
problems there is considerable literature which it is not our purpose to review 
here. 

1.3 The Burgers equation as testing ground for Navier— 
Stokes 

The Burgers equation, because of its known solutions, is frequently used for test- 
ing numerical schemes, particularly those intended for compressible flow (many 
of the Google hits are of this kind). If one is mostly interested in turbulence, 
as is the case for participants of the present School, Burgers equation turns out 
to be quite useful for testing - and mostly discarding - certain types of theories 
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of turbulence. Indeed, there have been many attempts to tackle the problem 
of the statistical theory of turbulence by adapting to it tools borrowed from 
field theory (for reviews, see Refs. [jD], [l8| [ll| pp[|). Such methods had little im- 
pact on the field until recently when they have permitted a real breakthrough 
in understanding the mechanism for intermittency and anomalous scaling (see, 
e.g. the lectures by G. Falkovich et al. in the same volume). In the past such 
field-theoretic methods have frequently involved formal expansions in powers 
of the nonlinearity, with Feynman graphs used for the bookkeeping of all the 
terms generated after averaging over Gaussian initial conditions and/or random 
forces. Since the Burgers equation has the same type of nonlinearity as the 
Navier-Stokes equation such methods are typically also applicable to the Burg- 
ers equation. Hence it is possible to find what they predict for the latter and 
to compare the results with those obtained by more reliable methods. From 
this point of view, that is of using the Burgers equation as testing ground, it 
is desirable to know the answers to questions similar to those generally asked 
for Navier-Stokes turbulence. For example, what are the scaling properties of 
structure functions; what are the probability distribution functions (pdf) of ve- 
locity increments and velocity gradients? Such questions will be at the center 
of these lectures. Whenever possible we shall comment on the corresponding 
Navier-Stoke issues. The emphasis will be exclusively on what happens in the 
real space-time domain in the limit of vanishing viscosity, which is of course not 
the same as naively putting the viscosity equal to zero. A number of interesting 
questions, requiring a finite viscosity, such as the pole decomposition pll 
will thus be left out. 



2 Basic tools 

In this section we introduce various analytical, geometrical and numerical tools 
which are useful for constructing solutions to the decaying (unforced) Burg- 
ers equation (|l|). Mostly, we shall deal with the deterministic equation, while 
making occasional comments on consequences for burgulence. 



2.1 The Hopf— Cole transformation and the maximum rep- 
resentation 

If in (0) with F = we set ip — 2v In 9 we obtain the d-dimensional heat equation 

d t e = v\i 2 e, (3) 

which can be solved explicitly if there are no boundaries. One thus obtains 



ij)(r,t) = 2i/ln 



I 1 

\ (47T^) d/2 



exp 



l r ~ a r 

2t 



d d a , (4) 



G 



where ?/>o( a ) is the initial potential. The limit of vanishing viscosity (y — ► 0), 
obtained by steepest descent, has the following "maximum representation" 



V>(r, t) = max hfo(a) ^ — J • (5) 

Note that the operation of taking a maximum is global in nature, whereas the 
viscous Burgers equation is a local partial differential equation. If ipo(a) is 
differentiable (i.e. the initial velocity uo (a) exists as an ordinary function rather 
than a distribution), the maximum in (g|) will be achieved at one or several 
points a where the gradient of the r.h.s. vanishes, that is, where 

r = a + tv (a) (6) 

In other words, r is the position at time t of the fluid particle starting at a 
and retaining its initial velocity vo(a). Hence, we can interpret a and r as be- 
ing, respectively, Lagrangian and Eulerian coordinates. Along this Lagrangian 
trajectory, the velocity being conserved, we have 

v(r,t)=v (a). (7) 

The map a^r defined by (|J) is called the naive Lagrangian map. It is not 
necessarily invertible: if there are several Lagrangian locations satisfying (|J) for 
a given r the only acceptable one is that which maximizes the argument on the 
r.h.s. of (^). As long as the Jacobian of the naive Lagrangian map @ 

J ( a '^= det (^'-^) (8) 

does not vanish the map is guaranteed to be invertible and the solution of the 
Burgers equation cannot have a singularity. For sufficiently smooth initial data 
with bounded second derivatives of V'o the first singularity appears at 

t* = TT7 vf ' ( 9 ) 

max a [A (a) J 

where A(a) is the largest eigenvalue of the Hessian matrix d 2 ipo/daida,j. 

In one dimension, we denote the velocity by u. Now, the time t* is the inverse 
of the absolute value of the most negative initial velocity derivative duo(a) / da. 
It is the first time at which the characteristics x — a + tuo(a) of the hyperbolic 
inviscid Burgers equation are crossing (Fig. 0). The first singularity in one or 
more dimension, is known as a "preshock" p3[ and plays an important role in 
the theory of pdf for velocity gradients and densities (Sections ^ and . 

Note that, for Gaussian random initial conditions, is itself random and 
can become arbitrarily small with very small but nonvanishing probability. As a 
consequence, most averaged quantities (e.g. the two-point correlation function) 
will have an essential singularity at t = 0. Note also that the distribution 
of eigenvalues of the Hessian matrix extends to infinite values in any finite 
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by initial velocities 



x 



Figure 3: Characteristics for the unforced one-dimensional Burgers equation in 
the (x, t) plane. 

dimension d, but becomes compactly supported (on a semi-circle) when dividing 
the eigenvalues by \fd and letting d — > oo. This is indeed a consequence of the 
properties of large random symmetric matrices, called Wigner matrices (see, 
e.g., Ref. H). 

2.2 Shocks in one dimension 

After the time the Lagrangian map @ ceases to be invertible. More precisely, 
for a given Eulerian position r there is more than one Lagrangian position a 
satisfying ([j|). This implies the presence of shocks in the Eulerian velocity 
field. In this subsection we consider the one-dimensional case and give various 
geometrical constructions of the solution (including shocks). 

First, let us define in the (x, V')-plane the Lagrangian manifold (a curve in 
one dimension) 

x = a + tuo(a) (10) 
ip = ipo(a) - ^ul(a), (11) 

where the second line is just the r.h.s. of (^]) without the maximum, evaluated at 
the (naive) Eulerian position a + tuo (a) . Fig. |] (upper) shows this Lagrangian 
manifold after the time t±. Hence, above some Eulerian locations x there is 
more than one branch and cusps are present at Eulerian locations such that the 
number of branches changes. Clearly, the correct Eulerian potential is obtained 
by taking the maximum, i.e. always the highest branch. Note that this potential 
will have one or several points with discontinuous slope, the right derivative 
being always greater than the left one. Hence the velocity, which is the negative 
space derivative of the potential (shown in the lower part of Fig. ||) will have 
discontinuities at shock locations with u- > u+. It is also possible to directly 
construct the velocity starting from the Lagrangian manifold in the (x, u)-plane 

x = a + tuo(a) (12) 
u = uo(a). (13) 
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If there is a single shock present, it follows obviously that its position is deter- 
mined by a Maxwell rule: the hashed loops shown on Fig. |^ (lower part) right 
and left of the shock should have equal areas. A Maxwell rule construction can 
become very cumbersome if there are several shocks present. 





Figure 4: Lagrangian manifolds for the potential in the (x, -0)-plane (upper) 
and the velocity in the (a;, it)-plane (lower). The heavy lines correspond to the 
correct Eulerian solutions. The vertical dashed lines delineate the multivalued 
region. 

Another geometrical construction uses the Lagrangian potential 

a 2 

<^(M)=# (a)-y, (14) 

whose negative gradient is obviously the naive Lagrangian map. We can rewrite 
(§ as 

x 2 

ttp(x, t) + — = max [tp(a, t) + ax] , (15) 

2 a 

which represents the potential as, basically, a Legendre transform of the La- 
grangian potential. (Note that the Legendre transformation is also used in the 
theory of multifractals.) The r.h.s. of ( [l5|) is equivalent to finding the largest 
algebraic vertical distance between the graph of the Lagrangian potential and 
the line of slope — x through the origin. If the graph is convex (second derivative 
negative everywhere), the maximum is attained at the unique point where the 
derivative has the value —x. Otherwise, it suffices to replace the graph of if by 
its convex hull <p c , that is the intersection of all half-planes containing the graph. 
This is illustrated in Fig. ^, which shows both regular points (Lagrangian points 
which have not fallen into a shock) and one shock interval, situated below the 
segment which is part of the convex hull. Again, it is possible to work directly 
with the (negative) derivative of the Lagrangian potential, namely, the naive 
Lagrangian map. The convex hull construction becomes then a Maxwell rule as 
shown in Fig. [| From this one can easily show that the speed of a shock is the 
half-sum of the velocities immediately to the right and to the left. 
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Figure 5: Convex hull construction in terms of the Lagrangian potential. 




Maxwell rule 



Figure 6: Naive Lagrangian map before (left) and after (right) appearance of a 
shock. The correct Lagrangian map is obtained by a Maxwell rule. 



Finally, the maximum formula (|5|) yields directly a "parabola construction" , 
illustrated in Fig. ^: a parabola with apex at x and radius of curvature pro- 
portional to t is moved down until it touches the initial potential V'o(a) at the 
Lagrangian location associated to x (or at two such locations if there is a shock). 
Which of the five geometrical methods given is more convenient depends on the 
application considered. The parabola construction is best for understanding 
evolution in time (cf Section |). It may be used, for example, to show that the 
long-time Eulerian solution has a sawtooth structure with shocks separated by 
ramps of slope 1/f (see Fig. |l^). The ramps are associated to high local maxima 
in the potential tpo- 

With random and homogeneous initial conditions there will be shocks (dis- 
continuities) at random Eulerian locations which do not cluster (unless we use 
non-smooth initial conditions as in Section From this it is easily inferred 
that, for p > 0, the structure functions 



S p (Ax, t) = (\u(x + Ax, t) - u(x, t)\ p ) 
behave, for small Ax, as 

S p (Ax,t)~C p \Ax\ p + C' p \Ax\, 



(16) 



(17) 
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Figure 7: Parabola construction of the solution. 



where the first term comes from regular (smooth) parts of the Eulerian velocity, 
while the second comes from the 0(|Aa;|) probability to have a shock somewhere 
in an interval of Eulerian length \Ax\. For < p < 1 the first term dominates 
as Aa; — > 0, while, for p > 1, it is the second. Hence, S p ~ | A.cc| ^ , with the 
exponents £ p as shown in Fig. |q. There are also higher-order corrections to 
the simple scaling law given in ( |17| ) which cannot be obtained by such simple 
arguments [jj0|. Note that a second-order structure functions with a behavior 
oc \Ax\ at small distances implies an energy spectrum E(k) oc k~ 2 as k — > aa. 



Figure 8: Exponent of the structure function of order p. Note the "phase 
transition" at p = 1. 

The "phase transition" at p = 1 seen in Fig. || is due to the isolated character 
of the dissipative structures (the shocks) , a feature not present in incompressible 
three-dimensional Navier-Stokes turbulence. 

2.3 Convex hull construction in more than one dimension 

Some of the methods used for the onc-dimensional case are readily extended to 
dimensions d > 1, for example the construction from the Lagrangian manifold 
in the (d + 1) dimensional space (x, tp). In Section ^ we shall use the multidi- 
mensional generalization of the convex hull construction, which we now briefly 
outline. We define the Lagrangian potential 



1 




1 



p 



(p(a,t) = — — +tip (a). 



(18) 
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and find, from (||), that 

Irl 2 

tip(r, t) + L±- = max [<p(a, t) + r • a . (19) 

2 a 

As before, this involves a (multidimensional) Legendre transformation which 
leads us to the construction of the convex hull in a (d + 1) dimensional space of 
the graph of the Lagrangian potential. In more than one dimension, singularities 
of convex hulls are considerably more involved. As a consequence, the equivalent 
of shocks are discontinuities across (d — 1) manifolds, but there are many other 
singularities of higher codimcnsion (the codimension is d minus the dimension 
of the object). 

In two dimensions the convex hull consists generically of four kinds of objects: 
(i) parts of the original graph, (ii) pieces of ruled surfaces, (iii) "kurtoparabolic 
points", to which we shall come back, and (iv) triangles (see Fig. ^|). The 
associated Eulerian objects are, respectively, (i) regular points, (ii) shock lines, 
(iii) end points of shocks and (iv) shock nodes. Likewise, in three dimensions we 
have two-dimensional shock surfaces meeting in triples at shock lines, meeting in 
quadruples at shock nodes. (Nodes are always connected to shock lines and never 
isolated.) Note that the Eulerian part of Fig. |^ looks just like a thermodynamic 
phase diagram, with the three shock lines playing the role of the liquid-gas, 
liquid-solid and solid-gas transition lines, the node playing the role of the triple 
point and the end point the role of the critical point. This is not accidental. 
In thermodynamics, equilibrium states are obtained by minimizing the Gibbs 
potential. This is equivalent to taking a Legendre transform of the internal 
energy in which the pressure and the temperature play the role of the Eulerian 
coordinates |^6| . This analogy holds also in higher dimensions: the classification 
of "Legendrian singularities" can be used both for studying the Burgers equation 
[ p7j and for studying multi- variable phase transitions [ ^8| . 

A more complete description of singularities is obtained by considering the 
metamorphoses of singularities as time elapses. A complete classification in 
two and three dimensions may be found in the appendix (supplement 2) by 
V.I. Arnold, Yu.M. Baryshnikov and LA. Bogaevski of Ref. ||. 

For random initial conditions the structure functions scale the same way as 
in one dimension. For example, the probability of having a (d— l)-dimensional 
shock intercepting an Eulerian segment of small length r is 0(r). The higher 
codimension structures give only subdominant corrections. 

2.4 Remarks on numerical methods 

Here, we give just some indications on how the (decaying) Burgers equation can 



be solved numerically. (For more details, see Refs. [J7j, |29|, 30 or any textbook 
on numerical methods for nonlinear hyperbolic equations; for the case with 
forcing, see Refs. |l5[ |3^] and references therein.) First, one can of course, 
solve the Burgers equation with viscosity. This should be avoided unless one 
is interested in what happens at dissipative scales (e.g., inside shocks). For 
the inviscid limit and only in the decaying case it is possible to construct the 
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Figure 9: Construction of the convex hull for a two-dimensional Lagrangian 
potential and associated Eulerian picture (Figure adapted from Ref. @). 



solution at time t directly from the initial condition without recourse to any 
time marching. One way is to directly use the maximum representation (||), 
assuming that Lagrangian and Eulerian locations have been discretized on the 
same grid. Then, for a given Eulerian r one searches the Lagrangian a which 
maximizes the r.h.s. If there are N grid points, this seems to require 0(N 2 ) 
operations, but it can actually be done in 0(N log 2 N) operations R ^9|. Such 
a strategy must be combined with suitable interpolations to increase accuracy 
and avoid getting complete garbage for derivatives p0| . In one dimension one 
can also use Lagrangian strategies with particle and/or shock tracking. To be 
consistent with the inviscid limit, the particles must stick upon collisions, for 
sufficiently smooth initial data, it may be possible to construct the solution 
from the Lagrangian manifold (|Io|)-(|ll|) or its multidimensional generalization, 
by just searching the maximum, for a given x, of the finitely many branches 
present. 

3 The Fourier— Lagrange representation and arte- 
facts 

In this section we show that formal manipulations of the inviscid Burgers equa- 
tion with random initial conditions, even though they include apparently terms 
of all orders, can nevertheless lead to completely incorrect results, e.g. for the 
energy spectrum. This section is entirely based on work by Fournier and Frisch 
[ p3[ . The theory is given in one dimension but similar results can be established 
in higher dimensions. 

In one dimension, it follows from that the Eulerian solution to 
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the initial value problem for the decaying Burgers has the following implicit 
representation: 

u{x, t) = uo(a) 

x = a + tua(a). (20) 

This becomes explicit if, instead of working with u(x,t), we use its spatial 
Fourier transform (27r-periodicity is assumed for convenience) 

1 f 27r 

u(k,t) = — e~ lkx u{x,t)dx (21) 
2^ Jo 

and make the change of variables x i— * a, to obtain 

1 f 2rr Bt 
u(k,t) = — e- ikx{a ^u (a)—da, x(a,t) = a + tu (a). (22) 
2tt J da 

Eq. (^2|) is called the Fourier-Lagrangian representation. A first integration by 
parts yields 

Mk,t) = —\ [ " e- lk{a+tUo(a)) u'{a)da. (23) 
27T ik J 

A second integration by parts leads then to 

u(k,t) = f e ~ lk(a+tUo{a)) da, k^O. (24) 

2tt ikt J 

If we now take random homogeneous Gaussian initial conditions, we can easily 
calculate moments of u(k,t) because they just involve averages of exponentials 
having the Gaussian initial velocity in their arguments. For example, the energy 
spectrum, related to the correlation function by 

(u(k,t)u(k',t))=E(k,t)5 k , k ,, (25) 

where 5k,k' is a Kronecker delta, has the following expression 

£( M )=^_l r e -ikh e -i k H*s 2 (h,o) dK (26) 

27T kH z J 

where S2(h,0) = (\uo(h) — uo(0)] 2 ^ is the second-order structure function of 
the initial velocity field. If the latter is smooth, as we shall assume, we have 
S2(h,0) oc h 2 for h — > 0. It then follows by a simple Laplace-type asymptotic 
expansion of (|2^) that 

E(k,t)(xk- 3 when k -> oo. (27) 

This is obviously the wrong answer: for Gaussian initial conditions there will 
be shocks with a non- vanishing probability for any t > 0. Their signature is a 



fc law in the energy spectrum at high wavcnumbers, as shown in Section 2.2 
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What went wrong? After the appearance of the first shock the Lagrangian 
map a <— > x is not monotonic and the change of variable from (|2l]) to ( p2| ) is 
valid only outside of the Lagrangian shock interval. Hence, in Q22| ) we should 
excise this interval from the domain of integration. If we do not remove it, we are 
actually calculating the Fourier transform of a function obtained by superposing 
the threee branches shown in Fig. [l^ with a plus sign for the two direct branches 
and a minus sign for the retrograde branch (the sign comes from the lack of an 
absolute value on the Jacobian dx/da in (E2|)). Obviously, this superposition 
has two square-root cusps as shown in Fig. This produces fc -3 / 2 tails in the 
Fourier transform and, hence, explains the spurious fc~ 3 energy spectrum. Note 
also that this superposition of three branches is not a solution to the Burgers 
equation, the latter being nonlinear. This phenomenon is not related to the 
well known non-uniqueness of the solution to the Burgers equation with zero 
viscosity without proper additional conditions [24] . 



Figure 10: Spurious solution of Burgers equation when three branches of a 
multi- valued solution are combined into one. 

The problem is actually worse than suggested so far. It is easily shown 
that if the the initial velocity is deterministic and smooth, the function of the 
time defined by (p2|), for fixed wavenumber k, is entire, that is, its Taylor series 
around t = has an infinite radius of convergence. There is no way to see the 
time i* of the first preshock from this function. A preshock is indeed an "ultra- 
violet" singularity which is not seen in the temporal behavior of a single spatial 
Fourier component. This result has an important consequence for the case of 
random Gaussian initial conditions. Suppose we simply ignore the viscosity in 
the Burgers equation and expand the solution to all orders in a temporal Taylor 
series around t = and then calculate various correlation functions and use 
Feynman graphs for bookkeeping of all the terms generated from averaging. We 
then find that the whole set can be resummed exactly and gives a spectrum 
with a fc~ 3 tail. Of course, the origin of the "resummation miracle" is the 
Fourier-Lagrangian representation. 

4 The law of energy decay 

An important issue in burgulence and turbulence is that of the law of decay 
at long times when the viscosity is very small. Before turning to the Burgers 
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equation let us recall a few things about the Navier-Stokes case. It is generally 
believed that high-Reynolds number turbulence has universal and non-trivial 
small-scale properties. In contrast, large scales, important for practical applica- 
tions such as transport of heat or pollutants, are believed to be non-universal. 
This is however so only for the toy model of turbulence maintained by prescribed 
large-scale random forces. Very high-Reynolds number turbulence, decaying 
away from its production source, and far from boundaries can relax under its 
internal nonlinear dynamics to a (self-similarly evolving) state with universal 
and non-trivial statistical properties at all scales. Karman and Howarth plf , 
investigating the decay of high- Reynolds number, homogeneous isotropic three- 
dimensional turbulence, proposed a self-preservation (self-similarity) ansatz for 
the spatial correlation function of the velocity: the correlation function keeps 
a fixed functional shape; the integral scale L(t), characteristic of the energy- 
carrying eddies, grows in time and the mean kinetic energy E(t) = u 2 (t) decays, 
both following power laws; there are two exponents which can be related by the 
condition that the energy dissipation per unit mass |i£(£)| should be propor- 
tional to u 3 /L. But an additional relation is needed to actually determine the 
exponents. The invariance in time of the energy spectrum at low wavenumbers, 
known as the "permanence of large eddies" (ll], [34| can be used to derive 
the law of self-similar decay when the initial spectrum E (k) cx k n at small 
wavenumbers k, with n below a critical value equal to 3 or 4, the actual value 
being disputed because of the "Gurbatov phenomenon" (see the end of this sec- 
tion). One then obtains a law of decay E(t) cx £- 2 (™+ 1 )/( 3 +' 1 ). (Kolmogorov 
[ |32| proposed a law of energy decay u 2 (t) cx i~ 10 / 7 , which corresponds to n = 4 
and used in its derivation the so-called "Loitsyansky invariant" , a quantity ac- 
tually not conserved, as shown by Proudman and Reid Q.) When the initial 
energy spectrum at low wavenumbers goes to zero too quickly, the permanence 
of large eddies cannot be used, because the energy gets backscattered to low 
wavenumbers by nonlinear interactions. For Navier-Stokes turbulence the true 
law of decay is then known only within the framework of closure theories (sec, 
e.g., Ref. |0|). 

For one-dimensional burgulcnce, many of these questions are completely 
settled. First, we observe that the problem of decay is quite simple if a finite 
spatial periodicity is assumed. Indeed, eventually, all the shocks produced will 
merge into a single shock per period, as shown in Fig. The position of the 
shock is random and the two ramps have slope 1 /t, as is easily shown using the 
parabola construction of Section |7|. Hence, the law of decay is simply E{t) cx 
t~ 2 . Nontrivial laws of decay are obtained if the burgulence is homogeneous in 
an unbounded domain and has the "mixing" property (which means, roughly, 
that correlations are decreasing with separation) . The number of shocks is then 
typically infinite but their density per unit length is finite and decreases in 
time because shocks are constantly merging. The E(t) cx i- 2 (™+ 1 )/( 3 +™) l aw 
mentioned above can be derived for burgulence from the permanence of large 
eddies when n < 1 iQ. For n = 0, this i~ 2 / 3 law was actually derived by 
Burgers himself fl35| . 

The hardest problem is again when permanence of large eddies does not 
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Figure 11: Snapshot of solution of decaying burgulence at long times when 
spatial periodicity is assumed. 



determine the outcome, namely for n > 1. This problem was solved by Kida 
[|| (see also Refs. §, || |§). 

We now give some key ideas regarding the derivation of Kida's law of energy 
decay. We assume Gaussian, homogeneous smooth initial conditions, such that 
the potential is homogeneous. Since a homogeneous function is not, in general, 
the derivative of another homogeneous function, we assume that the initial 
energy spectrum 

E Q {k)cxk n , n>l; fc^O. (28) 

This condition implies that the mean square initial potential J k~ 2 Eo(k) dk has 
no infrared (small-fc) divergence (the absence of an ultraviolet divergence is 
guaranteed by the assumed smoothness). 

A very useful property of decaying burgulence, with no known counterpart 
for Navier-Stokes turbulence, is the relation 

£(*) = J^>, (29) 

which follows by taking the mean of (0) in the absence of a driving force. Hence, 
the law of energy decay can be obtained from the law for the mean potential. 
The latter can be derived from the cumulative probability of the potential which, 
by homogeneity, does not depend on the position. By (JbJ), its expression at x — 
is 



Prob {Potential < tp} = Prob jVa, ip (a) < — + tpj . (30) 

Expressed in words, we want to find the probability that the initial potential 
does not cross the parabola a 2 /(2t) + tp (see Fig. |l2] ). Since, at long times t, 
the relevant tp is going to be large, the problem becomes that of not crossing 
a parabola with small curvature and very high apex. Such crossings, more 
precisely the upcrossings, are spatially quite rare. As a consequence of the 
mixing property, for long t, they form a Poisson process |57| for which 

Prob {no crossing }~ e -< N ®\ (31) 

where (N(t)) is the mean number of upcrossings. By the Rice formula (a con- 
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Figure 12: An initial potential which is everywhere below the parabola a 2 /(2t) + 
The probability of such events gives the cumulative probability to have a 
potential at time t less than 



sequence of the identity 5(Xx) = (l/|A|)<5(a;) 
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where H is the Heaviside function and 



m(a) =i/3 (a)-—. 



(32) 



(33) 



Since if>a(a) is Gaussian, the r.h.s. of (32) can be easily expressed in terms of inte- 
grals over the probability densities of ipo(a) and of d?po(a)/da (as a consequence 
of homogeneity these variables are uncorrelated and, hence, independent). The 
resulting integral can then be expanded by Laplace's method for large t, yielding 



(N(t)) ~ tV-ty 



t — * oo. 



(34) 



When this expression is used in (|31j) and the result is differentiated with respect 
to ip to obtain the pdf of p(ip), the latter is found to be concentrated around 
= (lni) 1 / 2 (see Fig. O). It then follows that, at large times, we have Kida's 
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Figure 13: A sketch of the pdf of the potential at long times, 
log-corrected \/t law for the energy decay 
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The Eulerian solution, at long times, has the ramp structure shown in Fig. [U] 
with shocks of typical strength u(t) = E 1 / 2 ^), separated typically by a distance 
L(t). The growth in time of L(t) takes place because correlated particles, which 
initially cannot be much apart, may propagate to far- apart locations at long 
times. 



u 










u(t) 






L(t) 









Figure 14: The Eulerian solution at long times t. The ramps have slope 1/t. In 
time- independent scales, the figure would be stretched horizontally and squeezed 
vertically by a factor proportional to t. 



The fact that Kida's law is valid for any n > 1, and not just for n > 2 as 
thought originally, gives rise to an interesting phenomenon now known as the 
"Gurbatov effect" : if 1 < n < 2 the long-time evolution of the energy spec- 
trum cannot be globally self-similar. Indeed, the permanence of large eddies, 
which is valid for any n < 2 dictates that the spectrum should preserve exactly 
its initial C n k n behavior at small wavenumbers k, with a constant-in-time C n . 
Global self-similarity would then imply a t~ 2 ( n + 1 ) / (s+n) j aw £ or ^ e ener gy de- 



cay, which would contradict Kida's law. Actually, as shown in Ref. 34 , for 



1 < n < 2 there are two characteristic wavenumbers with different time depen- 
dences, the integral wavenumber k^t) ~ (L(i)) -1 and a switching wavenumber 
k s {t) -C fcx,(t) below which holds the permanence of large eddies. It was shown 
that the same phenomenon is also present in the decay of a passive scalar |38| . 
Whether or not a similar phenomenon is present in three-dimensional Navicr 
Stokes incompressible turbulence or closure models thereto is a controversial 
matter |§, g|. 

For decaying burgulence, if we leave aside the Gurbatov phenomenon which 
does not affect energy-carrying scales, the following may be shown. If we rescale 
distances by a factor L(t) and velocities amplitudes by a factor u(t) = i? 1 / 2 (t) 
and then let t — ► oo, the spatial (single-time) statistical properties of the whole 
random velocity field become time-independent. In other words, there is a self- 
similar evolution at long times. Hence, dimensionless ratios such as the velocity 
flatness 

F{t) = (36) 

have a finite limit as t — oo. A similar property holds for the the decay of 
passive scalars [Q. We do not know if this property holds also for Navier- 
Stokes incompressible turbulence or if, say, the velocity flatness grows without 
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bound at long times. 



5 One-dimensional case with Brownian initial 
velocity 

Burgers equation, when the initial velocity is Gaussian with a power-law spec- 
trum oc k~ n , is what cosmologists call scale-free initial conditions (see Refs. 
H). Here, we consider the one-dimensional case with Brownian motion (in the 
space variable) as initial velocity, corresponding to n = 2. The general case, 
including higher dimensions, is discussed in Ref. Hi (an example of a 2-D sim- 
ulation with scale- free initial data is shown in Fig. |2|). 

Brownian motion is continuous but not differentiable (see Fig. |l5| ); hence, 
shocks appear after arbitrarily short times and are actually dense (see Fig. |l6| ) . 
Numerically supported conjectures made in Ref. || , have led to a proof by Sinai 
[|U of the following result: in Lagrangian coordinates, the regular points, that 
is fluid particles which have not yet fallen into shocks, form a fractal set of 
Hausdorff dimension 1/2. This implies that there is a Devil's staircase of di- 
mension 1/2 in the Lagrangian map (see Fig. |l8|). Note that when the initial 
velocity is Brownian, the Lagrangian potential has a second space derivative 
which is delta-correlated in space; this can be approximately pictured as a situ- 
ation where the Lagrangian potential has very strong oscillations in curvature. 
Hence, it is not surprising that very few points of its graph can belong to its 



convex hull (see Fig. 17) 



We will now give some highlights of Sinai's proof of this result. For this 
problem, it turns out that the Hausdorff dimension of the regular points (deter- 
mined in Ref. |Q) is also equal to its box-counting dimension, which is easier 
to determine. One obtains the latter by finding the probability that a small La- 
grangian interval of length £ contains at least one regular point which belongs 
simultaneously to the graph of the Lagrangian potential if and to its convex 
hull. In other words, one looks for points, such as R, with the property that the 



i. 




Figure 15: A realization of the Brownian motion curve. The parabola shows 
the root-mean-square velocity tx a 1 / 2 . 
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Eulerian velocity (h=0.5) 
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Figure 16: Snapshot of the velocity resulting from Brownian initial data. Notice 
the dense proliferation of shocks (from Ref. 0). 



(a) fP t (a) 



Figure 17: Sketch of the Lagrangian potential together with its convex hull 
(from Ref. @). 



graph of if lies below its tangent at R (see Fig. [19]). Sinai does this by the box 
construction with the following constraints on the graph: 

Left: the graph of the potential should be below the half line T_, 

Right: the graph of the potential should be below the half line T + , 

!1: enter (AF) with a slope larger than that of T_ by O^ 1 / 2 ) 
2: exit (CB) with a slope less than that of T + by 0(£ 1/2 ) 
3: cross (FC) and stay below (ED). 

It is obvious that such conditions ensure the existence of at least one regular 
point. (Move (ED) down parallel to itself until it touches the graph.) Note that 
A and the slope of (AB) are prescribed. Hence, one is calculating conditional 
probabilities; but it may be shown that the conditioning is not affecting the 
scaling dependence on £, 
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Figure 18: Left: the Lagrangian map looks like a devil's staircase. Right: 
standard devil's staircase over the triadic Cantor set, which is constant almost 
everywhere, except on a fractal (from Ref. Rj). 




Figure 19: The box construction used to find a regular point R (point of tan- 
gency with the graph entirely on one side of the tangent) within a Lagrangian 
interval of length £ (from Refs. || §). 

As the Brownian motion uo( a ) is a Markov process, the constraints Left, Box 
and Right are independent and hence, 

P reg (£) = Prob {regular point in interval of length £} 

= Prob {Left} x Prob {Box} x Prob {Right} (37) 

The scales of the box were chosen so that Prob {Box} is independent of I: 

Prob {Box} (38) 

Indeed, Brownian motion and its integral have scaling exponent 1/2 and 3/2, 
respectively and the problem with £ « 1 can be rescaled into that with I — 1 
without changing probabilities. 

It is clear by symmetry that Prob {Left} and Prob {Right} have the same 
scaling in i. Let us concentrate on Prob {Right}. We can write the equation for 
the half line T+ in the form 

F+ : a ^ V {a 2 ) + 5£ 3/2 + (d aV {a 2 ) + 7 ^ 1/2 ) (a - a 2 ), (39) 
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where 7 and S are positive 0(1) quantities. Hence, introducing a = a — the 
condition Right can be written to the leading order as 

.,2 



(u (a) + 7<? 1/2 ) da + Si 3 ' 2 + 



a 



> 0, for all a > 0. 



(40) 



By the change of variable a = j3i and use of the fact that the Brownian motion 
has scaling exponent 1/2, one can write the condition Right as 

(u (a)+^)da > -6, for all (3 G [O,^ 1 ]. (41) 



Without affecting the leading order, one can replace the Brownian motion by a 
stepwise constant random walk with jumps of ±1 at integer a's. The integral 
in (pi] ) has a geometric interpretation, as highlighted in Fig. 20, which shows a 
random walk starting from the ordinate 7 and the arches determined by succes- 
sive zero-passings. The areas of these arches are denoted S*, Si, —S n , S**. It is 
easily seen that 




Figure 20: The arches construction which uses the zero-passings of a random 
walk to estimate the integral of Brownian motion (from Refs. jl2| fjj). 



Prob{Rigiit} - Prob^i >0,S 1 + S 2 > 0, Si + ... + S n > 0} , 



(42) 



where n = 0(^ -1 / 2 ) is the number of zero-passings of the random walk in the 
interval [0, i^ 1 ]- The probability J42] ) can be evaluated by random walk methods 
(sec, e.g., Rcf. Q, Chap. 12, Sec. 7), yielding 



Prob {Right} - Prob{n first sums > 0} oc n~ 1/2 oc £ 1/4 . 



(43) 



By ([37]) , ( p8| ) and ( p3| ) , the probability to have a regular point in a small interval 
of length £ behaves as I 1 / 2 when i — * 0. Thus, the regular points have a box- 
counting dimension 1/2. 

6 Preshocks and the pdf of velocity gradients in 
one dimension 



In this section we shall determine the tail behavior of the probability density 
function (pdf) of the velocity gradient for one-dimensional decaying burgulence. 
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To explain some of the motivations for this study, it is useful to make a digression 
concerning the forced one-dimensional Burgers equation: 



d t u + ud x u = vd xx u + f(x, t), (44) 
u(x,t ) = u Q (x). (45) 

The latter displays much richer features than the unforced problem. The case 
where the force is random has often been studied as a prototyp e for a wide range 



of problems in non-equilibrium statistical physics (see Section 1.2). 

Eq. (^J) can also be used in the same spirit as the forced Navier-Stokes 
equation, namely to investigate universality of various statistical properties with 
respect to the forcing. For Navier-Stokes turbulence, when the force is confined 
to large spatial scales and the Reynolds number is very high, small-scale (inertial 
range) statistical scaling properties are generally conjectured not to depend 
on the forcing, except through overall numerical factors. Similar conjectures 
have been made for burgulence with large-scale forcing. For example, there is 
little doubt that, because of the presence of shocks, structure functions of order 
p > 1 have universal exponents equal to unity, as in the decaying case (see, e.g., 
Refs. jl(| Q). More controversial is the tail behavior of the probability density 
function (pdf ) of velocity gradients and velocity increments in the limit of zero 
viscosity when the force is a white-noise process in time. For increments, the 
problem was addressed for the first time by Chekhlov and Yakhot [^5), who 
considered a force with a power-law spectrum, acting both at large and at small 
scales. Concerning the pdf p(£ ) at large negative gradients £, it is generally 
believed that it follows a power law 

p(0«l?r s forC^-oo, (46) 

but the conjectured values of a differ markedly. Polyakov and Boldyrev [^7j , 
using a field-theoretical operator product expansion, predicted a = —5/2; E et 
at |44j], using a semi- heuristic approach in which preshocks (nascent shocks) 
are key, predicted a — —7/2; Gotoh and Kraichnan [Q, using a Fokker-Planck 
equation approach, predicted a — —3; more recent work by Kraichnan p9| 
favored a = —7/2. E and Vanden Eijnden |so[ developed a probabilistic 
formalism adapted to solutions with shocks and giving insight into many aspects 
of the problem; they proved that a < —3, and made a good case for a = 
—7/2. The question of the correct law for the case of white- noise forcing remains 
however open (we shall come back to this later). 

Actually, there is a situation much simpler than that originally considered 
in Ref. f44| , for which the argument in favor of a = —7/2 can be made rigorous, 
namely decaying burgulence. This closes our digression; in the remainder of this 
section we concentrate on the the unforced problem 

d t u + ud x u = vd x u, (47) 

in the limit of vanishing viscosity v — > and we follow Refs. |52|, [53). We 
assume a random initial velocity uq ~ —(dipy/da), deriving from a smooth 
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initial potential. Homogeneity is not required. The value a — —7/2 for the 
exponent of the pdf at large negative gradients is easily understood in this case. 
It is just the signature of the preshocks, the cubic root singularities in Eulerian 
coordinates, which appear when new shocks are created (see Fig. |T]). Preshocks 
constitute discrete events in space-time, contrary to shocks which persist in 
time (until they merge). These preshocks are the only structures giving large 

♦ 



u(x,t) 




X 



Figure 21: Eulerian structure of the solution just before a preshock, at the time 
ti, of a preshock and just after. 

finite negative gradients: shocks give infinite negative gradients (unless a finite 
viscosity is introduced) and the gradients in the immediate spatial neighborhood 
of a mature shock are not particularly large. A simplified presentation is given 
hereafter for the case of a single preshock; the contributions of several preshocks 
to the pdf are just additive. 

Let us suppose that the initial gradient dug /da has a minimum at a = 
(corresponding to an inflection point with negative derivative of the initial 
velocity), so that a shock will appear at time t = £* = —l/((duo/da)(0)) and 
at x = t*Uo(0). Without loss of generality, we assume uq(0) — (otherwise we 
perform a Galilean transformation to bring it to zero) . As the initial velocity is 
supposed to be sufficiently smooth, we can perform a Taylor expansion of the 
initial potential in the neighborhood of a — 0. We then have, locally, 

ipo(a) = c\a 2 - c 2 a 4 + h.o.t., (48) 

where c\ and c 2 are positive (random) constants and "h.o.t." stands for higher- 
order terms. The Lagrangian potential is locally 

a 2 t 
V = -y + #o(a) = -a 2 - tc 2 a A + h.o.t., (49) 

where r = (t — £*)/£*. The Lagrangian map outside the shock is thus 

x(a, t) = -d a ip(a, t) = -ra + 4tc 2 a 3 + h.o.t. (50) 

The Lagrangian potential, together with its convex hull, are shown in Fig. ^2|. 
It is convex for t < t*. At t = i*, there is a degenerate maximum with quartic 
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behavior, and, immediately after t* (for r > 0), convexity is lost and a shock in- 
terval is born. Given the symmetry, resulting from our choice of coordinates, the 
convex hull contains a horizontal segment extending between the two maxima 
a± = ±(t/(4c2)) 1 / 2 . The velocity gradient can be written locally as 



<P(a,t) 




T<0 

Figure 22: Normal form of the Lagrangian potential in the neighborhood of a 
preshock in one dimension. At the time of the preshock (r = (t — t*)/t* = 0), 
the Lagrangian potential changes from a single extremum to three extrema and 
develops a non-trivial convex hull (shown as a dashed line). 



dM*,t)=^M= y± 2 , (5i) 

O a X(a, t) -T + 12C2GT 

where a is the unique preimage of x by the naive Lagrangian map outside of the 
shock interval ]a_,a+[. Since, by (|50|), the relation between x and a is cubic 
at r = 0, the velocity gradient d x u(x,t ± ) cx |a;|~ 2 / 3 , which is unbounded. For 
any t ^ t+, the gradient remains bounded, except at the shock location. For 
r < 0, just before creation of the shock, the cubic relation between x and a still 
holds, except in a region of Lagrangian width of the order of t 1 / 2 , and hence of 
Eulerian width ~ t 3 / 2 , where the relation becomes linear to leading order. 

The question is now: what is the fraction of Eulerian space-time where 
d x u < £, with £ a large negative number? Because of the cubic root structure, 
x must be in a small interval of width ~ |£|~ 3 / 2 - The time must be sufficiently 
close to £* for this interval still to be in the region of validity of the cubic relation, 
that is, within ~ |x| 2 / 3 ~ Hence, the relevant space-time fraction or, in 

other words, the cumulative probability to have d x u < £ is proportional to 
|£|-5/2_ -phis gives a pdf cx |£|~ 7 / 2 a t large negative £'s. 

Actually, there is another contribution, also proportional to |£| -7 / 2 stem- 
ming from a small time interval r ~ |cc| 2 / 3 ~ |^| 1 just after t* when small- 
amplitude shocks are present which have not yet completely destroyed the cubic 
root structure (see Ref. |32|). Similar arguments can be used to show that there 
are power-law ranges with exponent —7/2 and +1 in the pdf of velocity incre- 
ments for decaying burgulence [52| . 

The preshock argument has first been introduced phenomenologically in 
Ref. E3] to predict pdf's of velocity gradients and increments for the case of 
white-noise in time forcing at large scales. In principle, in the presence of forcing, 
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spatio-temporal accumulations of preshocks, invalidating the —7/2 law, cannot 
be ruled out. Nevertheless, numerical evidence in favor of the —7/2 law has 
been recently obtained by one of us (JB), using particle tracking simulations 
with a shot-noise approximation to white noise. 



7 The pdf of density 

In cosmological applications of the adhesion model/Burgers equation, it is of 
special interest to analyze the behavior of the density of matter, since the large- 
scale structures may also be characterized as mass condensations. In Eulcrian 
coordinates, the mass density p satisfies the continuity equation 

%> + V- (pv) = 0. (52) 

The initial density is denoted by po{a). 

The question we intend to address here is the behavior, in the limit of van- 
ishing viscosity and at large p's, of the pdf of mass density p(p), when the 
initial velocity is random and smooth (and not necessarily homogeneous) . This 
problem was studied in Ref. Q, where it was shown that density pdf's have 
universal power-law tails with exponent —7/2 in any dimension. This behavior 
stems from singularities, other than shocks, whose nature is quite different in 
one and several dimensions. (Similar results can in principle be obtained for 
velocity gradients and increments which are, however, not scalars in more than 
one dimension.) 

In one dimension, the pdf of the mass density at large arguments is basically 
the same as the pdf of gradients at large negative arguments. Indeed, it is easy 
to show that, for any x not at a shock location, 

p(x,t) = po(a) (1 - td x u(x,t)) , (53) 

where a is the preimage of x by the Lagrangian map || . If now po is bounded 
from below and above (e.g., for uniform po), the result of the previous section 
implies that, for p — > oo, the pdf p(p) of the mass density satisfies a p -7 / 2 law, 
which is again the signature of preshocks. 

The key to studying this problem in more than one dimension is the geomet- 
ric construction of the solution via the convex hull of the Lagrangian potential 



(see Section 2.3). Conservation of mass (p2J) implies that the density is given at 



regular points by 

/ ,x Po( a ) ,r A s 
P(X '* )= J(M)' (M) 

where J is the Jacobian of the Lagrangian map. (The density is infinite within 
shocks.) Since the Jacobian is (up to a factor (— l) d ) equal to the Hessian of the 
Lagrangian potential (determinant of the matrix of second space derivatives) , 
it follows that large densities are typically obtained only near parabolic points 
(where the Hessian vanishes). However, arbitrarily close to a parabolic point 
there are generically hyperbolic points where the surface defined by tp crosses its 
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tangent (hyper)plane and which, therefore, do not belong to its convex hull. Yet, 
there exist in general exceptional "kurtoparabolic" points which are parabolic 
and belong to the boundary of the set of regular points (kurtos means convex 
in Greek). Near such points, arbitrarily large densities are obtained. In one 
dimension, the only kurtoparabolic points are the preshocks which are discrete 
space-time events in both Eulerian and Lagrangian coordinates. In two and 
more dimensions, kurtoparabolic points are also born at preshocks but live in 
general for a finite time; they reside on manifolds of spatial dimension (d — 1) 
(see Fig. 0). In Eulerian space, they are associated to boundaries of shocks (e.g. 
end points of shock lines for d = 2). 




Figure 23: Lagrangian potential in two dimensions with (a, b) coordinates, just 
after a preshock (left) and in the immediate neighborhood of a kurtoparabolic 
point (right). Continuous lines: separatrices between the regular part and the 
ruled surfaces of the convex hull; dotted-lines: vanishing of the Jacobian of the 
Lagrangian map. A and A' are a pair of kurtoparabolic points born with the 
shock. 



The determination of the large-p tail of the cumulative probability distri- 
bution of the density, P > (p), is equivalent to finding the fraction of Eulerian 
space-time where p exceeds a given value (see Ref. ]53| for details). The latter 
is determined by changing from Eulerian to Lagrangian coordinates and Taylor- 
expanding to the relevant order the Lagrangian potential near a kurtoparabolic 
point, in a suitable coordinate frame: 

^(a, t) ^Of + E [-f a) + foala 3 ] . (55) 



From (55), it is then easy to determine explicitly the line of vanishing Jacobian, 
the separatrix of the convex hull and the area where the density exceeds the 
value p (as illustrated in Fig. ^for the 2-D case). 

When p — ► oo, the cumulative probability can be estimated as follows 

P>(p) oc p~ 3/2 x p- 1 x 1 x . . . x I x ^ . (56) 

from a i from a 2 from a 3 ...a d from time 
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Figure 24: Projection in the two-dimensional case of the neighborhood of a 
kurtoparabolic point. 



Hence, the cumulative probability P > (p) oc p~ 5 / 2 in any dimension; so that the 
pdf of the mass density has a universal power-law behavior with exponent —7/2. 
We have seen that the theory is rather different in one dimension and higher 
dimensions, because kurtoparabolic points are persistent only in the latter case. 
However, the scaling law for the resulting pdf is the same in all dimensions. 
Actually, two orthogonal spatial directions, a\ and 0,2 in (|5^), play the same 
role as space and time in one dimension. 

It is now clear that, for burgulence, the algebraic tails of the pdf of ve- 
locity gradients or of the density stem from singularities. Turning briefly to 
incompressible three-dimensional Navier-Stokes turbulence, we note that mea- 
surements of pdf 's for space or time derivatives of Eulerian velocities have not 
yet revealed power-law tails, but such tails may just have been, so far, "lost 
in the experimental noise". There has indeed been considerable speculations 
about singularities of the Navier-Stokes equations in the inviscid limit pj| . If 
singularities with divergent gradients are present, they will give power-law tails, 
at least as intermediate asymptotics when the viscosity is small (the converse 
is however not true, since statistical effects not related to singularities can also 
give power laws) . The confirmed absence of power laws would probably rule out 
singularities. 



8 Kicked burgulence 

8.1 Forced Burgers equation and variational formulation 

In the limit of vanishing viscosity and when no force is applied, the Burgers 
equation just means that fluid particles keep their initial velocity until they 
stick together in a shock. So, until merger, the position X(t) of a given fluid 
particle will depend linearly on time: 

X(t) = X(t ) + (t-t )u (X(t )), u(X(t),t) = u (X(t )). (57) 
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Figure 25: Trajectories of particles in the decaying case (a), in the continuously 
forced case (b) and in the kicked case (c). The thick lines are shock trajectories. 



When a force is applied, fluid particles trajectories, before merger with a shock, 
follow forcing-dependent continuous trajectories governed by 



dt 2 



X(t) = f(X(t),t), -X(t) = u(X(t),t); 



(58) 



thus their dynamics can be rather complex (see Fig. g5|) . 

Following Ref. |30), we shall be concerned here with the one-dimensional 
case where the force is a sum of impulses (or kicks), concentrated at discrete 



times tj s: 



E 

3>3o 



f(x,t)=Y / f j (x)6(t-t j ), 



(59) 



where tj — to is the initial time, and tj +\ is the time of the first kick. The 
kicking times tj's and the kick fj(-)'s are prescribed. They can be either fixed 
or random. The meaning of such a forcing is that, between kicks, we let the 
solution evolve as a solution of the unforced problem. At each kicking time tj, 
we discontinuously change the velocity field by the amount fj(x): 



u(x, t 



3 + ) 



u(x,tj-) + fj(x). 



(60) 



This is an intermediate case between decay and time-continuous forcing. Such 
forcing implies a piecewise-linear time dependence of the position of a given 
fluid particle (see Fig. p5| ). 

It is of interest to notice that this kind of discrete-in-time forcing can be 
applied also to the Navier-Stokes equations, with features of decaying turbulence 
still present to some extent. The original motivation for introducing such a 
forcing was to approximate white-noise-in-time forcing by discrete random noise, 
also called shot noise. But actually, the kicked case displays interesting features 
of its own. As will be shown later, the problem can be understood in terms of 
area-preserving mappings to which we can appl y K AM theory (see Ref. p4| and 
references therein) and Aubry-Mather theory [|55[ p6[ . 

We will focus on the space-periodic case. Namely, we assume that both 
the initial condition uq(-) and the kicks fj(-) are periodic with period 1, with 
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respect to the space variable. For the moment, let us also assume that the initial 
velocity and the kicks both have zero spatial mean value over the space period 
[0, 1[. Since the mean velocity is conserved by Burgers dynamics, we have 



u(x, t)dx = 



(61) 



at all times. This constraint implies that the velocity potential tp(x,t), defined 
by u[x, t) = —d x ip(x,t) 7 is itself periodic in space. Let us define the kicking 
potentials Fj(-), so that 

f j ( x ) = -±F j (x). (62) 

It is then easy to write the potential at any time t, using between successive 
kicks the standard maximum representation (^) for decaying solution in the 
limit of vanishing viscosity (this is reexpressed here as a minimum in order to 
minimize a suitable action function), to obtain 



ip(x,t) = -: 



2(t-tj) 



Hyj,tj-) - Fj{yj) 



(63) 



where the index J is such that tj < t < tj + i. Repeating this step as often as 
necessary to work our way back to the initial time, we obtain 



ip(x, t) = - min [A (x, t; {yj}) - ipo(yj )] , 

\yjlio<i<J 

where A is an action which has to be minimized, 

x2 J-i 



A{x,t-{y 3 }) 



[x - yjf 
2(t - tj 



E 

3= jo 



2(tj+i - tj) 



F 3+i(yj+i) 



(64) 



(65) 



There is a similar representation for the case where the forcing f(x, t) — —d x F(x, t) 
is continuously applied, namely 



tp(x, t) = - min [A (x, t; y(-)) - ipo(y(t ))] 
y(-) 



(66) 



The minimum is now taken over continuous curves y(-) such that y(t) = x, the 
action being given by 



A(x,t;y{-)) = 



l(y( s )f-F(y( S ),s) 



ds. 



(67) 



This representation goes back to work by Oleinik |57j on general conservation 
laws. It can be derived as the continuous limit of the discrete formulation 
when letting the time between kicks tend to zero. Many features of the forced 
Burgers equation were obtained by E et al. ]58| |. As we will see, the key notions 
introduced by E et al., such as minimizers, global minimizer and main shock, 
are still valid in the case of discrete-in-time forcing. 
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First, we will introduce the notion of minimizing sequence (or minimizer). 
In terms of fluid particles trajectories, the minimum representation ( p4| ) just 
means that, to obtain the solution at time t and at some Eulerian location 
x, one has to look at all possible trajectories reaching x, and choose between 
them those which minimize the action. The sequence for which the minimum is 
achieved is, by definition, a minimizer. In general, there is only one minimizing 
trajectory arising at a given x. But for a countable set of ir-values, there are 
several minimizing trajectories. These correspond to particles coalescing in a 
shock. 

A minimizer can be explicitly characterized by requiring the vanishing of 
the derivatives, with respect to all the j/j's, of the argument of the minimum in 
(|64|). A minimizing sequence then has to verify the following Euler-Lagrange 
equations: 

Vj+i = Vj+fjiyj), (68) 

Vj+i = yj + (tj+i-tj)l v j + fj{yj)}> ( 69 ) 

where Vj = {y.j — yj_i)/(tj — is the velocity at the location yj just before 
the kick. These equations have to be supplemented by the following initial and 
final conditions: 

v io = My 'j ), ( 7 °) 

x = y.j + (t -tj)v J+ i. (71) 

Note that Vj+i = u(x,t). The Euler-Lagrange map is area-preserving. It is 
also explicitly invertible, so that for a given (x, v), one can reconstruct the past 
history of a particle, except if a shock sits at x. 



8.2 Periodic kicks 

From now, we will focus on a particular case of forcing which displays globally 
the same features as random forcing but is much easier to handle. Namely, 
following Ref. [^0|, we consider the case of time-periodic kicks: the kicking 
potential is the same at each kick, Fj(x) — G(x) for all j, and the time interval 
is constant, tj+i —tj = X, for convenience. The force can then be written 

f(x,t)=J2g(x)5(t-j), (72) 

3>j0 

where g = —dG/dx. 

We now show, following Ref. ||(| , that the solution to the Burgers equation 
with this kind of forcing converges exponentially fast in time to a periodic solu- 
tion Uoo {x,t). Snapshots of the time-periodic solution for one instance of kicking 
are shown in Fig. ^6|; Fig. ^ shows the exponential relaxation to u QO (x, t). 

Actually, the convergence to a unique solution at long times is related to 
properties near a fixed point of the two-dimensional dynamical system defined 
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Figure 26: Snapshots over one time period of the velocity for the limiting so- 
lution corresponding to the kicking force g(x) — s'mx — cos(3ic) on the space 
period [0, 2ir[ (see upper inset). The main shock is located around x = ir; the 
global minimizer, here a fixed point, is the point of vanishing velocity common 
to all curves. Notice that during each period, two new shocks are born and two 
mergers occur. 



by the Euler-Lagrange map which reads here 

Vj+i = Vj + g(yj), (73) 
Vj+i = Vj + Vj+g(yj). (74) 

A fixed point (y*, v*), obviously, satisfies = and = 0. The latter 

expresses that the kicking potential achieves an extremum at i = y±. Let 
P = (x c ,0) be the particular fixed point of the map (|7^)-([74[), which corre- 
sponds to the location where the forcing potential achieves its maximum over 
the space period. This point is hyperbolic because the linearized system in its 
neighborhood has two real eigenvalues A > 1 and 1/A, where 



i d 



■2 



A = 1 + c + Vc 2 + 2c, c = --_G(a: c ). (75) 

2 dx z 

In the phase space (x,v), two globally invariant curves are associated to the 
corresponding eigendirections. These are (i) the stable manifold r^ s ', associated 
to 1/A, which is the set of points (x,v) converging to the fixed point under 
iteration of the map (because the eigenvalue is less than one), and (ii) the 
unstable manifold T^, associated to A, and generated by inverse iteration (see 
Fig. ^8|). An arbitrary continuous curve in the (x,v) plane which intersects the 
stable manifold will, under iteration, converge exponentially fast to the unstable 
manifold at the rate 1/A. 

In the language of Burgers dynamics, the curve in the [x, v) plane defined by 
an initial condition uq(x) will be mapped after some kicks into a curve very close 
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2 4 6 8 10 12 14 

number of kicks 



Figure 27: Exponential relaxation to a time-periodic solution for the same 
forcing as in Fig. |2(| with three different initial conditions, as labeled. 
J Q v \u(x, n_) — Uoo(x, l_)|dx/(27r) is plotted vs the number of kicks. 




Figure 28: Sketch of a hyperbolic fixed point P with its stable (r( s )) and its un- 
stable (r<")) manifolds. A curve, which intersects T^ s \ will eventually converge 
to r^ 1 - 1 under iteration of the map. 

to the unstable manifold. To understand this mechanism of convergence, let us 
take an initial time to tending to — oo and look at the behavior of the solution 
at time t = 0. The trajectory of the hyperbolic fixed point P corresponds 
to the so-called global minimizer. The global minimizer is the trajectory of 
a fluid particle never to be absorbed by a shock. Such a global minimizer is 
unique, and every minimizing trajectory converges exponentially fast to the 
global minimizer as t — > — oo |5q| . This is illustrated in Fig. [29|a. By definition 
of the unstable manifold, each point (jjj,Vj) of a minimizer belongs thus to 
r(") and every regular part of the graph of the limiting velocity belongs to the 
unstable manifold. Now we turn to the construction of the main shock. Since 
A > 0, every minimizing trajectory starting from a point on the right (resp. 
left) of the global minimizer approaches it as t — ► — oo from the right (resp. 
left). Hence, there exists x r (resp. x{), the rightmost (resp. leftmost) location 
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(a) 




(b) 



Figure 29: (a) Minimizers on the (x,t) cylinder; initial time to is taken at — oo. 
Shock locations are characterized by having two minimizers (an instance is at 
xi); the fat line X — X q IS the global minimizer. (b) Unstable manifold I" 1 -"' on 
the (x,v) cylinder which passes through the fixed point P = (x c ,0); the bold 
line is the graph of the limiting periodic solution. The main shock is located at 
xi = x r , and another shock at x\ corresponds to a local zig-zag of between 
A and B. 



from which a minimizer approaches the global minimizer from the right (resp. 
left). By periodicity in space and uniqueness of the global minimizer, these two 
points are actually the same: x r = xi mod 1. If we shift the periodicity interval 
to [jcjjav], we can draw on the (x, u)-cylinder. The regular parts of the 
limiting solution belong to this graph. By construction there is thus a shock at 
x = xi — x r (see Fig. p9[b). This is the main shock, the unique shock which 
exists for an infinite time. In Burgers dynamics, shocks are born and then they 
may merge. The main shock is a shock which has always existed when letting 
the initial time tend to — oo. The other shocks are associated to the regions 
where r^™' is multi-valuated in x. Their locations are determined by requiring 
that the action be the same at points such as A and B in Fig. [2S|b. 

8.3 Connections with Aubry— Mather theory 

So far, we have exclusively considered zero-mean-value initial conditions. Let 
us briefly consider the case where 

/ u(x,t)dx — / uo(x)dx = a > 0. (76) 
Jo Jo 

The Burgers problem is then in exact correspondence with the description of 
equilibrium states of the Frenkel-Kontorova model [^|. In the latter, one has a 
one-dimensional chain of atoms connected by elastic springs in the presence of 
a space-periodic potential. The potential energy, which must be minimized to 
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obtain the (classical) ground state, is 



#({%}) 



G(x), 



(77) 



where a is the unstrctched distance between atoms. This problem was inves- 
tigated by Aubry 55 and Mather fjjfj. The representation (77) matches the 
action minimizing representation for Burgers equation with a mean velocity a. 
The connection between the forced Burgers equation and Aubry-Mather theory 
was investigated by Jauslin et al. |p0||, E and Sobolevski |62fl. 



0.2- 
o.l- 



f(x) = sin (2k x) 
< u > = 0.340492 




Fixed point with accumulation 
of shocks from the left 



0.2 0.3 0.4 0.5 



Figure 30: Velocity profile at the first transition to a rational nonvanishing 
rotation number when increasing the mean velocity a — (u) from 0. Note the 
accumulation of shocks. 

For a = 0, the global minimizer is a trivial ground state, associated to a 
fixed point, but for a ^ 0, it is much more complex. Within some intervals of 
the parameter a, the global minimizer lives on a periodic orbit associated to a 
rational rotation number p (asymptotic slope of the trajectory when t — > -co). 
The graph of p as a function of a is actually a Devil's staircase. The transitions 
between the intervals of the mean velocity corresponding to rational rotation 
numbers display interesting phenomena, such as accumulations of shocks (see 
Fig. 1). 
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